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^ ' Abstract 

ON ; 

We present a general method for extracting the Gell Mann-Low logarithmic derivative 
Ph! of an effective charge of an observable directly from data as a mean for empirically 

verifying the universal terms of the QCD /3-function. Our method avoids the bi- 
ases implicit in fitting to QCD-motivated forms as well as the interpolation errors 
introduced by constructing derivatives from discrete data. We also derive relations 
between moments of effective charges as new tests of perturbative QCD. 
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1 Introduction 

An effective charge M encodes the entire perturbative correction of a QCD observable; 
for example, the ratio of e + e~7* — > hadrons annihilation to muon pair cross sections 
can be written 

R ^_ w , <**£. - M — > = K ^_ w U + MVf) ) , (1) 

cr(e + e — ► /U + // ) y 7r y 

where R° e + e - is the prediction at Born level. More generally, the effective charge (Xa{Q) 
is defined as the entire QCD radiative contribution to an observable Oa{Q) 0: 

0,(A) = 0\ ( Sa + «±£>) , (2) 

where 5a is the zeroth order QCD prediction (i.e., the parton model), and au(A)/7r 
is the entire QCD correction. Note that 5a = or 1 depending on whether the 
observable A exists at zeroth order. Important examples with 8a — 1 are the e + e~ 
annihilation cross-section ratio and the r lepton's hadronic decay ratio, 

R = r(r~ -> t/ T + hadrons) = ^ / + a T (m r ) \ 
r(r~ -> z/ r e~z/ e ) r y % J ' 

In contrast, the effective charge cty(Q) defined from the static heavy quark potential 
and the effective charge ct >2 jets defined from e + e~ annihilation into more than two 
jets, cr >2 jets ' nave $A — 0. 

One can define effective charges for virtually any quantity calculable in perturba- 
tive QCD; e.g., moments of structure functions, ratios of form factors, jet observables, 
and the effective potential between massive quarks. In the case of decay constants of 
the Z or the r, the mass of the decaying system serves as the physical scale in the 
effective charge. In the case of multi-scale observables, such as the two-jet fraction 
in e + e~ annihilation, the arguments of the effective coupling a2j e t{s, y) correspond to 
the overall available energy and characteristic kinematical jet mass fraction. Effective 
charges are defined in terms of observables and, as such, are renormalization-scheme 
and renormalization-scale independent. 

The scale Q which enters a given effective charge corresponds to its physical 
momentum scale. The total derivative of each effective charge (Xa{Q) with respect to 
the logarithm of its physical scale is given by the Gell Mann-Low function: 



* A [a A (Q, m , Q m] = ^' \ 4 

dlogQ 

where the functional dependence of ^a is specific to the effective charge a a- Here m 

refers to the quark's pole mass. The pole mass is universal in that it does not depend 

on the choice of effective charge. It should be emphasized that the Gell Mann-Low \^ 

function is a property of a physical quantity, and it is thus independent of conventions 

such as the renormalization procedure and the choice of renormalization scale. 

A central feature of quantum chromodynamics is asymptotic freedom; i.e., the 
monotonic decrease of the QCD coupling otA{Q 2 ) at large spacelike scales. The em- 
pirical test of asymptotic freedom is the verification of the negative sign of the Gell 
Mann-Low function at large momentum transfer, a feature which must in fact be true 
for any effective charge. 

In perturbation theory, 

^=-4°>^-4 i} 4-^ 2} 4+--- ( 5 ) 

At large scales Q 2 » m 2 , where the quarks can be treated as massless, the first 
two terms are universal || and basically given by the first two terms of the usual 
QCD j3 function for N c = 3 

4 1} = ^ = -+-iv« (6) 

A 8 4 12 F ' A K ' 

Unlike the /5-function which controls the renormalization scale dependence of bare 
couplings such as aj^ s (fi), the ip function is analytic in Q 2 /m 2 . In the case of the ay 
scheme, the effective charge defined from the heavy quark potential, the functional 
dependence of N F y(Q 2 /m 2 ) is known to two loops ||. 

The purpose of this paper is to develop an accurate method for extracting the 
Gell Mann-Low function from measurements of an effective charge in a manner which 
avoids the biases and uncertainties present either in a standard fit or in numerical 
differentiation of the data. We will show that one can indeed obtain strong constraints 
on ^ A and ^ A from generalized moments of the measured quantities which define 
the effective charge. We find that the weight function /(£) which defines the effective 
charge ^/(A) from an integral of the effective charge cxa(Q) can be chosen to produce 



maximum sensitivity to the Gell-Mann Low function. As an example we will apply 
the method to the e + e~ annihilation into more than two jets. Clearly one could 
also extract the Gell Mann-Low function directly from a fit to the data, but the 
fact that we are dealing with a logarithmic derivative introduces large uncertainties 
0. Our results minimize some of these uncertainties. In addition, our analysis 
provides a new class of commensurate relations between observables which are devoid 
of renormalization scheme and scale artifacts. 

One can define generalized effective charges from moments of the observables. 
The classic example is a T (A) where A is the generalization of the lepton mass. The 
relevant point is that R T can be written as an integral of R e + e ~ f|, as follows: 



i? T (A 5 



£/tf 



A2 ds 



S(i-iV(i+£ 



A 2 



A 2 J 



A 2 



Rp 



(«), 



(7) 



where q/ are the quark charges. As a consequence of the mean value theorem, the 
associated effective charges are related by a scale shift 

a T (A) = a R (yfi = A T ), (8) 



The ratio of scales A T /A in principle is predicted by QCD ||: The prediction at NLO 
is [HI 



A 



cxp 



19 __ 169 a R (A 7 
24 
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(9) 
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Such relations between observables are called commensurate scale relations (CSR) |J. 
The relation between R T and R e + e - suggests that we can obtain additional useful 
effective charges by changing the functional weight appearing in the integrand. Indeed 
it has been shown |7j that, starting from any given observable O a we can obtain new 
effective charges a A f by constructing the following quantity 



O ai (A) = cJ a 



%;:Mt)^ 



(10) 



where C is a constant and /(£) is a positive arbitrary integrable function. In order 
for Oaj to define an effective charge a A f through 



Af (A) 



o° Af [s A 



a Af (AY 



TX 



(11) 



it is necessary that Ai(A) = A]A and A2(A) = A2A, with both Ai and A2 constant. 
Then, by the mean value theorem, a A f is related again to a A by a scale shift 



«a/(A) = a A (A Af ), 



(12) 



with Ai < A A j < A2. An important observation [0] is that PQCD predicts A a/ = 
Aaj/A to leading twist. If we ignore quark masses so that the two first coefficients of 
the Gell-Mann Low function are constant, one has 



a A (>/s) a A (A) 



TX 

1 
+ 4 



TX 



*> 2 [ T . 



#0, ( s 

T ln U 



'a A (AY 



+ 



TX 



(13) 



2^1n(- 



'a A (A)' 



TX 



If we now use eqs.([TUp and flTT|), we find |7| 



a Af {A) a A {A) 



TX 



TX 



2 I 0f 



'a A {A)' 



TX 



+ - 



,2 J 2/ 



^-^-2*i-^ 



1/ 



'0/ 



'0/ 



'au(A)' 



TX 



(14) 



where Iif = J x i /(£)(ln£ 2 ) d£ 2 is independent of the choices of observable A and scale 
A, but only provided that A X (A) = AiA and A 2 (A) = A 2 A. Replacing s by A\^ in eq. 
( |T3"| ) and comparing with eq. ([T4]), we find 



A 



a/ = exp 



! i/ 



^0 

2/ 0/ 4 



ill 






a A (A) 



TX 



(15) 



In general the commensurate scale relation will have the following expansion 

(n) (0< A (A)\ n 



ln\ Af (A) = Y,a) n 



ra=0 



TX 



(16) 



where the first three coefficients are independent of A. Note that the above formulae 
are only valid inside regions of constant Np and sufficiently apart from quark thresh- 
olds. If we include the mass dependence, the effective charges, by the mean value 
theorem, are still related by a scale shift, although it cannot be written in the simple 
form of eq. (|i~5|) . Indeed, even the lowest order of X A f would have a small dependence 
on the energy and the effective number of flavors appearing in \^o- 



2 Obtaining the Gell Mann-Low function directly 
from observables 

The main practical obstacle in determining the Gell Mann-Low function from exper- 
iment is that it is a logarithmic derivative. One can try to obtain the value of the 
parameters of the \1/ function from a direct fit to the data using the QCD forms, but 
any approximation to the derivative of the experimental results implicitly requires 
extrapolation or interpolation of the data. In order to observe a significant variation 
of the effective charge a a one needs to compare two vastly separated scales. This is 
illustrated in Fig.l. However, to approximate \P(\/5) — Aa / i(v / s)/(Aln y/s) with a 
huge separation between y/s and y/s' is not very accurate since then the value for 
Aa^/A In y/s is the slope of the Q straight line in Fig. 1 instead of that of P, which 
gives an O (A In y/s) 2 error. If we want to obtain $ from a finite difference approxi- 
mation, we need to interpolate A In y/s — > 0, but in this case the experimental errors 
will most likely be much larger than the required precision. Such an interpolation 
procedure has already been applied in ref. [Q near the r region to test the running of 
a s (including appropriate corrections to the leading twist formalism). In this energy 
region the value of the QCD coupling is rather large, and the interpolation yields 
evidence for some running. However, it has also been pointed out in ||, that the 
value of the coupling extrapolated from the r region to high energies appears small 
compared to direct determinations. 

In the next section we shall use the effective charge formalism to derive several 
expressions within leading twist QCD which relate the intrinsic \I/ a function of a^ 
directly to the observables O a- We shall show that with just three data points we 
can obtain good sensitivity to the value of \&o without any numerical differentiation 
or fit. 

2.0.1 Differential Commensurate Scale Relations 



Let us formally differentiate eq. ( |10T ) with respect to A 
dO Af (A) 1C 



d\ A 



Xif(X2)0 A (A 2 )-\if(X 1 )0 A (A 1 



20 Af (A) C f (W ds /^vg ^CvVA) ^ (17) 



A Ai ( A lA )^ A 2 " vv ' A d(yfs/K) 




Figure 1: Finite difference approximation of 'J. If one takes Aln-^/s very small, the errors can 
be larger than Aa, and the result will be meaningless. This can be avoided by choosing very far 
separated points y/s and y?, but then the approximation yields the slope of line Q instead of that 
of P. 

The first term in the right-hand side can be obtained directly from the data on O a . 
This is also the case for the second term, after using eqs. (Q) and fll2j) , since 



Af {A) = Af U A + 



«a/(A)~ 



Af [5 A + 



a A {A f Y 



§fO A (A,), 



18) 



7T I J \ IT 

Note that O a * and O a are known constants. Finally, there is a choice of /(£) which 
allows us to recast the third term in the right-hand side of eq. ( |TT| ) and provide a 
direct relation between the data and the effective charge. Namely, we choose 

df(0 



Z- 



d£ 



Pf(0, 



(19) 



with p any real number. That is, up to an irrelevant multiplicative constant, we take 

/(o = e- (20) 

With this choice eq. (|TT|) can be simply written as 



dO Ap (A) 2C 
dA A 



X p 2 +2 O a (A 2 ) - X^ 2 O a (A 1 ) 



P + 2 
A 



Ap (A). 



(21) 



Note that, to simplify the notation, we have substituted the / subscript by p. In 
terms of ty A this means 

da Ap (A) _ ttA dO Ap (A) 



*a p (A) = A 



dA 



O 



Ap 



dA 



(22) 



But using its definition, we can easily see that 

2CO° A 



O 



(p + 2; 



(Af-Af 2 ), 



(23) 



so that, using eq. (J18D , we arrive at 



*4p(A) 



71- 



p + 2 
O . 



,P+2 



P+2, 



A^O A (A 2 ) - X^OaJAi 
\p+ 2 _ y+ 2 



O a (A p ) 



(24) 



Note that we have just written ty Ap (A) directly in terms of observables. Therefore, 
we have related the universal \l/o an d ^i coefficients directly to observables, without 
any dependence on the renormalization scheme or scale. 

Up to this point A x and A 2 are arbitrary. In order to illustrate the meaning of 
eq.([24T), we now choose Ai = and A2 = 1, so that eq. (|24"D becomes: 



^Ap(A) 



n ^ [ ° a{A) " ° A{Kp) 



(25) 



Let us remark that, although it may look similar, the above equation is not the 
finite difference approximation 



MA) ^(A)-a,(A-AA) +0(AAi 



(26) 



which is a good numerical approximation to ^ A (A) when AA is very small. In 

contrast, eq. (|24]) , is exact (at leading twist) no matter whether A — A p is big or small. 

However, we do not want to set Ai = 0, since then the integrated effective charges 



defined in eq.(10), contain higher twist contributions which are unsuppressed at low 
energies, and our leading twist formulae would be invalid in practice. In addition, 
some observables like the number of jets produced in e + e~ annihilation are only well 
defined above some energy, which becomes a lower cutoff in the integral of eq. (|TDD . 

Nevertheless, by choosing A and A2 appropriately, we can obtain any value of 
Ai 7^ and A2 7^ 0, even if we set Ai = 1, and so we will do so in the following. That 
is: 



^Ap(A) 



7T- 



p + 2 

o A 



\ P 2 +2 Q A (A 2 ) - Q A (A) 
K +2 - 1 



O a {A, 



(27) 



which is an exact formula relating ^ A with the observable O a at three scales A < 



A p < A 2 



It happens, however, that we are interested in measuring not the ^Ap intrinsic 
function but *£?a itself. We thus arrive at our final result: 



V A (A\ P (A)) 



K 

1+ i 



n ■ 



p + 2 

o A 



X p 2 +z O A (A 2 ) - O a (A) 



O a {A\ 



{21 



Af 2 - 1 

where we have also defined X p = A p /A. Note that ^ a appears in the above equation 
both at A p and A through the A' coefficient, defined as dX/dLogA, which only van- 
ishes at leading order. Therefore, if we include higher order contributions the above 
equation is not enough to determine \I/ a at one given scale. 

Let us work out first the implications of eq.(^) at leading order, since it contains 
all the relevant features of our approach. 



2.1 Leading order 

Suppose then that we had three experimental data points at s a < s& < s c . In order 



to apply eq. fl27D , we first identify A 2 = Js c /s a and then we obtain the p such that 

^ = yW V 

The Ik p integrals are given by 



fc! 



l kp 



P/2 



7^ 

1 .7=1 



^^^(InA^H-iy 



s v(p/2+ m*- 07 (p/2+iy 



(29) 



Thus, at leading order we have to obtain p from 



ln^ = 21nA p = ^ 



,p/2+l 



ln(s c /s a ) 



' h P s p J 2+1 - s p J 2+1 P/2 + 1' 
which can be evaluated numerically. 

As we have already commented, at leading order A' = 0, and therefore 



(30) 



* 



Sb = 7T 



p + 2 

0°A 



SPJ 2+1 A {^) ~ sPj 2+1 A {^-q) 

o/2+l p/2+1 

oc i>a 



- Ca(v^) 



(31) 



Let us remark once more that these are leading-twist formulae, and s a , s&, s c should 
lie in a range where higher twist effects are negligible. 

2.2 Beyond leading order 

As we have already seen, if we go beyond the leading order contributions, we have to 
use eq. (]28|) , which does not completely determine the value of ^/a at a single scale. 



9 



In principle, we need an additional equation. In fact, the A' term can be neglected. 
Intuitively, this is due to the very slow evolution of a a- Let us give some numerical 



values; first, we will write 



K 



tf A (A)fy,(A), 



with 



fip(A)^ dlnXp 



2E 



na 



(n) 



' Oa(A) 
0° 



n-\ 



8, 



(32) 



(33) 



d(a A (A)) n= i \ ^ A 

From PQCD we know that the expansion of ^a starts with a\. Thus, the A' term in 
eq.(^8f) is an 0(a A ) effect. It should only be taken into account if we are interested 
in \l/ up to that order. Numerically, the expected value of ^(A) at the energies we 
will be using, ranges from 10~ 2 to 2 x 10~ 2 at most. In addition, Q ranges from 
3 x 10~ 2 to 0.5. Thus, even in the worst case, the A' term contribution would be 
slightly smaller than 1% of \l/. If that term is to be kept, then we need and additional 
equation involving a fourth data point. We have found that the final error estimate 
increases since it is much harder to accommodate four points sufficiently separated 
within a given energy range. It seems that 1% accuracy is the lower limit for this 
method. If additional higher twist corrections are included, it could be possible to 
extend the energy range to separate the points and improve the precision. 

Therefore, in what follows we will use eq. (J3ID . However, the NLO p parameter 
is now obtained by solving numerically the equation 



In* 



1 



„P/2+l 
o c 



-P/2+1 



+ 



^0 

2 



Sa P/2 + 1 

s a s c y/ 2+1 \n 2 (s b /s a ) 



(34) 



1 



r P/2+1 _ s p/2 + 1 )2 



'a 



0° A 



(P/2 + I) 2 . 

where s a < Sb < s c and ^/sb = \ p y/s^ and ^J!^ C = A2 \/s^- Note that now ^0 is an 
input, but the output is the NLO \1/ function. 



3 Error estimates 

Although they have inspired our approach, observables with 5a ^ are not well 
suited for our method, because the relative error in Ga(E) becomes at least one 
order of magnitude larger for the effective charge (Xa(E). For example, using the 
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e + e _ hadronic ratio defined in Sect.l, if we introduce a 1% error in R e + e -, the error 
in a.R is 0(20%) and we have to separate the data points over five orders of magnitude 
to obtain \& r with a 10% precision. In practice, that renders the method useless. 

The problem we have described is avoided if we use an observable with 5a = 0. 
That is the case, for instance, of the e + e~ annihilation in more than two jets, 
cr >2-iets( s '^) = °tot — diets' wnere V is used to define when two partons are unre- 
solved ||| (i.e. their invariant mass squared is less than ys). This process does not 
occur in the parton model since it requires, at least, one gluon. Note that ^ and ^ 
are independent of y. 

At LO we can work with exact results, but as soon as we introduce higher orders, 
there is some degree of truncation in the formulae. We have therefore first constructed 
simulated data following a model that corresponds to the exact LO equations. Let 
us remark that these are models, not QCD. They are obtained by the truncation of 
a a at a given order. Thus, in principle, they will have some different features from 
QCD, as for instance, some residual scale dependence. In the real world this will 
not occur. However, we have worked out these examples for illustrative purposes to 
obtain a rough estimate of the errors. 

3.0.1 Leading order 

What we call the LO model is to use 

a A (Q) _ a(M z ) V ln (Q 2 ) / « A (M Z ) \ 2 

7T '" 7T 2 VMf/ V 7T J ' l ' 

exactly. We have taken aA{Mz) as the reference value for simplicity. Note, however, 
that the derivative of the above expression is 

*„ = -£(«^)\ (36) 

which is a constant which differs by 0(a/n) s terms from the LO PQCD result 

.,(«) = -*»(«^))'. (37) 

In Table 1 we can see the estimates of the relative errors in our determination 
of ^a, which depend on the different position of the data points, as well as in their 
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v^ (GeV) 


v^ (GeV) 


v^ (GeV) 


ao a /o a 


A^a/^a 


30 


100 


300 


1% 


3% 


3% 


9.1% 


400 


640 


1000 


1% 


6.2% 


3% 


18.6% 


500 


875 


1000 


1% 


4.9% 


3% 


14.6 % 



Table 1: Estimated relative errors in the determination of ^o using the LO equations. We assume 
the relative error AOa/Oa in the measurements of Oa ■ The estimates correspond to an observable 
with a vanishing parton model contribution (5a = 0) such as e + e~ annihilation into more than two 
jets, cr >2 _jets- 

errors AOa- Since the observable vanishes in the parton model, the relative error in 
a a is exactly that of Oa- 

The results in the table deserve some comments. 



First, the values of y^ and y^ have to be chosen to maximize their distance, 
within a region of constant Np. Thinking in terms of 0" >2 _iets' they correspond 
either to the region where both energies are sufficiently above the b-quark pair 
threshold but still below ti production, or both are above the ti pair threshold, 
in regions accessible at NLC. 



• Second, we have chosen the same relative error for the measurements at the 
three points. The intermediate energy ^/sl is then tuned to minimize the error, 
which is obtained assuming the three Oa measurements are independent. 

Let us remark once again that we have not used at any moment the value of ^O) 
which is obtained from the data using this method. If we want to use higher order 
contributions, using the value of ^o as an input, we would obtain information about 
higher order coefficients, like \l/i if we were to work at NLO. 
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3.0.2 Beyond leading order 

The NLO model is now given by: 



a A (Q) a{M z ) * (Q 2 \(a A {M z ) 

■In 



1 

+ 4 

and therefore, we obtain 






,2t„2 



Z, 



'a A (M 2 



7T 



(38) 



MQH 4(-M) 2 _l!(-M) 3 , (39) 

which is the QCD NLO ^/^ result up to 0(a/7r) 4 terms. In contrast with the LO 
case, obtaining p now requires some truncation of the formulae when passing from 
eqs. ( |I3"D and ([H]) to eq. (|Hj). This is very interesting since we can thus obtain an 



estimate of the theoretical error due to truncation, which will be present in the real 
case too. It can be seen in Table 2 in the rows where AOa = 0, and it is usually 
0(1%). 

Again we have also considered the experimental AOa{E{) uncertainties. The final 
error given in the last column is estimated assuming that the four experimental errors 
and the one due to truncation are all independent. Note that when passing from a 1% 
experimental error to a 3%, the total error is not multiplied by 3, since the truncation 
error does not scale. 

The fact that we obtain larger errors in the NLO case may seem surprising, but 
it is not. The reason is that the LO is a very crude approximation of the \& a QCD 
scaling behavior. In the LO model, the \1/ function was a constant, but in the NLO it 
changes with the energy scale, as it occurs in the realistic case. Indeed, the evolution 
of a a at high energies becomes much slower so that the difference between a a at two 
given points is smaller at NLO than at LO. Hence, for the same relative errors, the 
relative uncertainties in the NLO \1/ function are much bigger. Of course, we expect 
the real data to show a behavior much closer to the NLO model. 

3.1 Using more than three points 

The advantage of fitting the data is that we can reduce the errors by larger statistics. 
But that is also true for our method. Up to now we have only used three points of 
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v^ (GeV) 


v^ (GeV) 


v^ (GeV) 


ao a /o a 


A^a/^a 


30 


100 


300 


0% 


2% 


1% 


2.7% 


3% 


7.5% 


400 


640 


1000 


0% 


.9% 


1% 


10% 


3% 


29% 


500 


875 


1000 


0% 


1% 


1% 


10% 


3% 


30% 



Table 2: Error estimates at NLO. 



data, but in the realistic case we expect to have several points at each energy range. 
It is then possible to form many triplets of data points, one at low energies (^/s^), 
another at intermediate energies (y/sb), and a last one in the highest range (\/sb)- 
Each one of these triplets will yield different values and errors for \1/, which can later 
be treated statistically, thus decreasing the error estimates given in Table 2. 



4 Conclusions 

We have obtained an exact and very simple relation between the Gell Mann-Low ^ 
function of an effective charge of an observable and its integrals. These results are 
renormalization-scheme and renormalization-scale independent. By choosing specific 
weight functions, these relations can provide an experimental determination of the 
PQCD \1/ function, thus testing the theory and setting bounds on the properties of 
new particles that would modify the expected QCD behavior. 

We have shown that a good candidate for this study is the e + e~ annihilation to 
more than two jets, since it is a pure QCD process. Even within the simple leading- 
twist formalism, which limits the applicability range, we have found that with just 
three precise measurements in present or presently planned accelerators, it could be 
possible to determine the \1/ function without making a QCD fit or any interpolation 
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and numerical differentiation of the data, eliminating the specific uncertainties of 
these methods. Thus we can obtain a determination of \& with different systematics. 
It also seems possible to extend the method and ideas, to include higher twist effects 
which will allow the use of a wider range of energies. This could result in an even 
more powerful set of tests of perturbative QCD. 
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